Polar bosons in one-dimensional disordered optical lattices 
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We analyze the effects of disorder and quasi-disorder on the ground-state properties of ultra-cold polar bosons 
in optical lattices. We show that the interplay between disorder and inter-site interactions leads to rich phase 
diagrams. A uniform disorder leads to a Haldane-insulator phase with finite parity order, whereas the density- 
wave phase becomes a Bose-glass at very weak disorder. For quasi-disorder, the Haldane insulator connects 
with a gapped generalized incommesurate density wave without an intermediate critical region. 



Introduction. The interplay between disorder and interac- 
tions plays a crucial role in the physics of strongly-correlated 
systems [1]. Disorder in non-interacting systems leads to An- 
derson localization |2||, which in one dimension occurs for 
vanishingly small disorder yfl. For the particular case of 
bosons in a lattice potential, interactions have been shown, 
both in ID 14J47D and higher dimensions lalSI], to induce in 
the presence of disorder a phase diagram characterized by 
three phases: a superfiuid (SF) phase, a gapless localized in- 
compressible phase known as Bose-Glass (BG), and a Mott- 
insulator (MI) occuring at commesurate lattice fillings. 

Ultra-cold atoms in optical lattices offer an extraordinarily 
controllable scenario for the detailed analysis of the compe- 
tition between disorder and interactions. Disorder in the on- 
site energies may be implemented in various ways in these 
systems, including the use of speckle II OH 1311 . binary mix- 
tures 111 4 - 11711 . and bichromatic combinations of two mutu- 
ally incommensurate lattices 111 811 . Recently, localization has 
been experimentally observed in non-interacting cold gases in 
ID and 3D speckle 119142111 . and bichromatic potentials 12211 . 
Bichromatic lattices constitute a peculiar type of disorder, 
rather a quasi-disorder, realizing the so-called Aubry-Andre 
model 1230 . The effects of interactions in ID lattices with 
quasi-disorder have been recently studied II24H27I1 . Particu- 
larly interesting is the existence of a gapped localized phase, 
the so-called incommensurate density wave (ICDW), which 
results from the quasi-periodicity of the potential. 

Polar gases are attracting a growing attention mostly mo- 
tivated by expe riments on atoms with large magnetic mo- 
ments 128143011 , and especially by recent groundbreaking 
experiments on polar molecules Bill . Due to the dipole- 
dipole interaction, these gases present an exceedingly rich 
physics I32U33II . Polar lattice gases are particularly interest- 



ing, mostly due to the qualitatively new features introduced by 
dipole-induced inter-site interactions 0411 . In particular, inter- 
site interactions may allow for the realization of the so-called 
Haldane-insulator (HI) phase 13511 . a gapped phase character- 
ized by a nonlocal string order parameter. 

In this Letter we show that the interplay between on-site 
and inter-site interactions and disorder leads to a rich physics 
for lattice bosons with nearest neighbor interactions. In par- 



ticular, in the presence of uniform disorder the HI phase is 
preserved, although with finite parity, up to a finite disorder, 
where a phase transition into a BG is produced. On the con- 
trary, in the presence of quasi-disorder the HI is connected, 
without any intermediate critical region, to a gapped general- 
ized ICDW phase occurring for non-polar gases. Other phases 
are discussed in detail for both types of disorder. 

Model. As mentioned above, the main qualitatively new 
feature of polar lattice gases concerns the significant inter-site 
interactions. Polar interactions between sites placed j > 
sites apart decay as 1/j 3 . Although interactions for j > 1 do 
play a role in the physics of polar gases, especially in what 
concerns the existence of crystalline phases at any fractional 
filling (Devil's staircase 06113711 ) for very weak on-site inter- 
actions and sufficiently large dipoles, the most relevant prop- 
erties of polar lattices gases may be understood from a model 
with only nearest-neighbor interactions: 
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with b\, bi the creation/annihilation operators for bosons at 
site i, tii = b}bi, t the hopping amplitude, U (V) the on- 
site (inter-site) interaction, and e, the on-site energy. 

In the following we consider an average filling factor 
n = l. Assuming small deviations from this average, we 
may introduce an effective spin via the standard Holstein- 
Primakoff (HP) transformation Sf = Srii = 1 — Ui, S^ = 
y/2 — riibi. The resulting spin Hamiltonian resembles to a 
large extent an antiferromagnetic spin-1 XXZ model with uni- 
axial single-ion anisotropy (the mapping is however imper- 
fect, especially at low U, due to occupations n > 2, and extra 
terms appearing when introducing the HP transformation): 



H = jJ2[sfs, 
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with J = 2t, D = U/U, and A = V/2t. For D > and 
A > three ground-state phases are possible 113 811 . a Neel 
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FIG. 1: Phase diagrams of bosons with nearest-neighbor interactions: (a) unperturbed case; (b) staggered on-site energy; (c) uniform disorder; 
(d) quasi-disorder. Figures (b)-(d) are obtained for U/t = 5. See text for details. 



antiferromagnet (AF) for sufficiently large A > 0, a large-Z? 
phase for dominant single-ion anisotropy, and the gapped Hal- 
dane phase 13911 characterized by a non-local string-order 114011 . 
Unraveling the HP transformation, the equivalent bosonic 
phases are: a density wave (DW) phase, characterized by 



are polarized by adding one more particle or eliminated by 
coupling two extra hard-core bosons at the chain edges in or- 
der to form a singlet state [46]. Both procedures destroy the 
energy gap for the MI and ICDW phases. We hence introduce 
two different definitions for the gap to the first excited state, 



a finite Odw = lim|i_ ]7 -|_ Hxl ((—l)* -J ' <5ni<Snj), and a static E G and E' G , which are respectively calculated without and 
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cently measured in site-resolved experiments |42j); and a 
Haldane-insulator (HI) phase, with finite string-order Os 
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renormalization group (DMRG H43H ) calculations for soft- 
core polar bosons 051 kill have shown that indeed a HI oc- 
curs in a window V C {1) (U) < V < V C (2) (U) ft At large 
U a direct first-order MI-DW transition occurs. All the pre- 
vious phases are insulating, being characterized by an expo- 
nentially decaying single-particle correlation function G(r) = 

(bj&i+r) / ' V '( n i) ( n i+r) fall . At low U occupations n > 2 
break down the boson-spin mapping, and a superfluid (SF) 
phase may be found, which is a Luttinger-liquid characterized 
by an algebraically decaying G(r) ex r~ 1 / 2K . A MI to SF 
Kosterlitz-Thouless (KT) transition occurs at sufficiently low 
U and V when the Luttinger parameter K > 2 l4lll45ll ). The 
MI-HI Gaussian transition line is a Luttinger-liquid with Lut- 
tinger parameter 2 > K C (U) > 1/2 yj]. At low U but large 
V our DMRG calculations show a DW to SF transition with a 
critical K = 1/2 ft (FigHJt). Below the HI region we find 
as well a SF region for K > 2. 

Numerical method. We employ below DMRG calcula- 
tions with open boundary conditions for determining the 
ground-state phases for different types of disorder. Special 
care must be taken with the edge states in the HI phase, which 



with the mentioned polarization/elimination procedure. For 
the MI and ICDW, E G ^ and E' G = 0; for HI and DW by 
Eg = and E' G ^ 0; and for the gapless phases (SF and BG) 
Eg = E' G = 0. The Luttinger exponent K of the SF phases 
is extracted by fitting the single-particle correlation function 
G(r) using finite-size conformal corrections IJ471 . The num- 
ber of particles is conserved and the maximal bosonic occupa- 
tion per site is n max = 8 for U ^ (14811 - We consider lattice 
sizes of TV = 55, 89, 144 and 233 sites (see discussion be- 
low). The number of optimal states kept in the DMRG ranges 
from 300 to 500. Statistical deviations are particularly rele- 
vant in the case of uniform disorder discussed below. For that 
case we have performed up to 60 realizations per case (we 
refer to the supplementary material for a discussion on the 
error bars). We have checked that our results converge to 
the exact values in the non-interacting limit (where we keep 
n m ax = N), and to earlier calculations for interacting sys- 
tems with V = Jg] (see the supplementary material). 

Staggered potential. As a starting point we first consider 
the case of a commesurate staggered on-site energy e, = 
( — 1) J A, which resembles the case of staggered magnetic 



fields in spin chains. For the latter case it is known (|49 



that the Haldane phase is continuously connected to the AF 
phase. Due to the explicit symmetry breaking, the AF phase 
is a singlet and not a doublet, and a Gaussian transition oc- 
curs between large-D and AF phases. We hence expect that 
the HI phase continuously connect to the DW phase for even 




FIG. 2: (a) Gap E' G as a function of A for the uniform disorder case 
with U — 5£, V — 3.3£, and different lattice sizes; (b) For the same 
parameters, string-order Os and parity-order Op, for N = 144. 




FIG. 3: Excitation gap of the ICDW as a function of A for the quasi- 
disordered case with U = 5£ and V — 0, and different lattice sizes; 
(inset) Static structure factor S(q) for A = 7£ (N = 144); S(q) 
presents side-peaks corresponding to the 1/(1-/3) quasi-periodicity. 

vanishingly small A, as observed in our numerical simula- 
tions (Fig. 01b)). Note that a finite SF region opens for the 
bosonic case between the MI and the DW phases. 

Uniform disorder. We consider next a uniform disorder, 
in which e, has a uniform random distribution in the inter- 
val [—A, A]. The presence of disorder induces a gapless lo- 
calized Bose-Glass (BG) phase 0, H, H], separating the MI 
and the SF phase. At unit filling, and intermediate U values, 
for growing A a MI-BG transition is followed by a BG-SF 
transition, and a final SF-BG transition at larger A , the SF- 
BG boundary forming a characteristic finger- like shape |]6| |7|] . 
The SF-BG boundary is characterized by a Luttinger parame- 
ter K = 3/2 J4J,|5[]. Our numerics recover these well-known 
results for V = (see Fig. Q2c) and the supplementary ma- 
terial). For growing V the intermediate SF region shrinks 
and eventually vanishes, whereas the MI phase stretches up 
to a critical V = V} , which marks the MI-HI transition for 
A = 0. Close to Vc numerical calculations cannot resolve 
whether a BG opens at the MI-HI boundary for vanishingly 
small disorder. Although RG calculations to the lowest order 
predict a stable MI-HI boundary at weak disorder 115 ill . RG cal- 
culations using Akp terms to the density 11521 15311 show that a 
BG should appear at the MI-HI boundary at vanishingly small 
disorder for K < 3/2 and at finite disorder for K > 3/2. 
For the value U/t — 5 in the figure, K = 1.18 and we hence 
expect a BG region all the way to vanishing disorder. 

we obtain that a gapped phase sur- 




For V c (1) < V < V c {2) 



FIG. 4: (a) Gaps Eg and E' G (extrapolated to infinite-size systems 
from our results for N = 55, 89 and 144) as a function of A for the 
quasi-disordered case with U — 5£, and V = 3.3£; (b) For the same 
parameters, string-order Os and parity order Op, for N = 144 



vives up to a critical A given approximately by the Haldane- 
like gap at A = (Figs.QJc) and[2|a)). This phase presents a 
finite Os- The survival of this gapped phase for finite uniform 
disorder resembles the persistence of the Halda ne p hase for 
weak random magnetic fields in spin-1 chains 15 ll l52l 15411 . 
This phase however is characterized by a finite Op which 
grows with increasing A (Fig. 0b)). Whereas in the pure 
HI phase Op = due to the fluid-like spatial distribution 
of defects 6n — ±1, the finite Op relates to the growing lo- 
calization of defects with increasing disorder, while still keep- 
ing string order. Resorting to the language of "solid-on-solid" 
models 14011 . we may map defects Sn = ±1 into spin-1/2 par- 
ticles, and sites with Sn = as empty sites. In this analogy the 
pure HI phase would be an AF-spin-ordered fluid, whereas the 
DW phase would be an AF-spin-ordered solid. The Haldane- 
like phase in the presence of disorder could be understood as 
an AF-spin-ordered glass-like in this analogy. 

On the contrary the DW is destroyed even for very weak 
disorder, similar to the destruction of the AF phase in spin- 1 
chains 15 ll 15211 . This results from the proliferation of domain 
walls in accordance with the Imry-Ma argument 115 511 , Our 

(2) 

results show for V > Vc a large reduction of the gap for 
A <C V, as well as the expected domain formation. However, 
the convergence of the gap towards zero with the system size 
is very slow, since the length of the Imry-Ma domains scales 
as l/\/A, and hence much larger lattice sizes are demanded 
to observe the gap destruction at very weak disorder. 

Quasi-disorder. We consider at this point the case of 
quasi-disordered potentials. Such potentials may be easily 
generated in ultra-cold gases by means of bichromatic lat- 
tices 122U56I1 . in which a second lattice incommensurate with 
the first lattice is added. In that case the system experiences 
an on-site energy tj — A cos(2ir/3j + <f>), in which j3 charac- 
terizes the incommesurability of both lattices, and A is given 
by the potential depth of the second lattice. In the following 
we consider fi = (V5 — l)/2. The chain length is chosen 
following the Fibonacci series N = 55, 89, 144, 233 in order 
that the potential is closest to be periodic with N. We have 
checked that in most cases the results are not affected by the 
choice of 0, and we hence set cf> = in the following. 



As mentioned above, it has been recently shown 124 12511 
that quasi-periodic potentials allow for an ICDW, which is a 
localized but gapped phase, contrary to the BG. The ICDW 
may be understood from the beating of both lattice frequen- 
cies which leads to super-wells with a characteristic length 
scale 1/(1-/3). A filling 1 — j3 (or /3) hence leads to the 
formation of an approximate density wave. For filling factor 
n = 1, a possibly gapped localized phase (i.e. a generalized 
ICDW) may occur as well, and a MI-SF-ICDW-SF-BG sce- 
nario cannot be excluded 02511 . Our results for V — show 
a clearly gapped localized ICDW whose gap converges to a 
finite value when extrapolating to infinite size (Fig. [3). This 
phase is characterized by a peak in the static structure factor 
S(k) given by the approximate 1/(1 — f3) periodicity (inset of 
Fig.0. Hence, at V = we observe a MI-SF-ICDW-SF-BG 
phase diagram (Fig.QJd)). The last SF region is rather narrow, 
and within numerical accuracy it is compatible with a critical 
line. At the SF-MI transitions we obtain numerically K ~: 1, 
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compatible with results obtained for spinless fermions 1157 1. 
For the SF-ICDW transitions we find as well K ~ 1. For 
V < Vc l) the MI-SF-ICDW-SF-BG topology is maintained, 
but the MI and low SF regions shrink vanishing at V c ■ 

For Vc < V < Vc and low A, we obtain as for the 
uniform disorder case a HI phase with growing Op when A 
increases (Fig. E][58]). This is again due to the growing lo- 
calization of defects Sn = ±1 by the quasi-disorder while 
keeping string-order. However, contrary to the uniform dis- 
order case, the quasi-disordered case allows for the general- 
ized ICDW phase which has also finite Os and Op. As a 
consequence, the increasing pinning of defects for increasing 
A leads to a connection between HI (E' G > 0) and gener- 
alized ICDW (Eq > 0), without any gapless region in be- 
tween (Fig.|4]i. Our DMRG calculations, including the deter- 
mination of fidelity susceptibility 15911 . show that the general- 
ized ICDW and the HI phases form an overall gapped region. 

(2) 

Finally, for V > V c we observe the opening of a finite DW 
region, which contrary to the uniform disorder case is not im- 
mediately destroyed at small A. The DW region undergoes at 
finite A an Ising transition into the generalized ICDW phase. 

Conclusions. In summary, uniform disorder and quasi- 
disorder lead to rather different phase diagrams for polar 
bosons in optical lattices. For uniform disorder a HI phase 
with finite parity survives up to a finite disorder where a HI- 
BG transition occurs, whereas the DW phase becomes a BG 
for any finite disorder. On the contrary, in the presence of 
quasi-disorder the HI phase continuously connects with a gen- 
eralized ICDW phase without any intermediate critical region, 
and the DW survives up to a finite disorder. The predicted 
phases may be detected using state of the art techniques. In 
particular, the string order may be detected using similar site- 
resolved techniques as those recently employed for the detec- 
tion of non-local parity order in MI phases 14211 . 
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In this supplementary material we briefly comment on the determination of the error bars for the case of uniform 
box disorder. 

We determine the SF-BG boundary by monitoring the behavior of the single-particle correlation function G(r), 
which has a power-law decay in the SF phase. This decay can be fitted very well using conformal field corrections 
for open boundary conditions, allowing for the calculation of the corresponding Luttinger parameter K . The SF-BG 
phase boundary is determined at the point in which K = 3/2 |l|. In order to benchmark our calculations we have 
established the phase diagram for disorded non-polar gases. Figure [1] shows that the calculated SF-BG boundary is 
in excellent quantitative agreement with that obtained in Ref. [2j using quantum Monte Carlo techniques. 

In order to determine the error bars, we have studied for different V and A up to 60 disorder realizations. We 
randomly pick different sub-sets of these realizations and average for each sub-set the correlation G(r). The error 
bars indicate for a given V the regime of A values in which for some, but not all, sub-sets correlations are best fitted 
by a power-law with K = 3/2. 

We establish the error bars between gapped and gapless phases, i.e. MI-BG and HI-BG using a similar procedure. 
We evaluate again a large number of disorder realizations for different lattice sizes. We perform finite-size scaling to 
obtain the excitation gap. The error-bar region is that at which some disorder realizations present finite gap whereas 
other realizations are not gapped. 
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FIG. 1: Phase diagram for the non-polar Bose-Hubbard model with uniform box disorder, for a chain length L — 89, with 30 
disorder realization per point. 



